Photon statistics dispersion in excitonic composites 



o 
o 

(N 

>v 



(N 



> 
On 
00 ■ 

in : 
o . 

o ■ 

9 L,: 
^— > ■ 

G ■ 

s : 



G. Ya. Slepyan and S. A. MaksimenkcQ 

Institute for Nuclear Problems, Belarus State University, 11 Bobruiskaya Str., 220050 Minsk, BELARUS 

Linear media are predicted to exist whose relative permiability is an operator in the space of 
quantum states of light. Such media are characterized by a photon statistics-dependent refractive 
index. This indicates a new type of optical dispersion - the photon statistics dispersion. Interaction 
of quantum light with such media modifies the photon number distribution and, in particular, the 
degree of coherence of light. An excitonic composite - a collection of noninteracting quantum dots 
- is considered as a realization of the medium with the photon statistics dispersion. Expressions 
are derived for generalized plane waves in an excitonic composite and input-output relations for a 
planar layer of the material. Transformation rules for different photon initial states are analyzed. 
Utilization of the photon statistics dispersion in potential quantum-optical devices is discussed. 

PACS numbers: 42.50.Nn, 42.50.Ct, 78.67.Hc, 42.50.Pq 



I. INTRODUCTION 

The search and investigation of novel physical pro- 
cesses transforming the photon number distribution is 
one of the most important problems of quantum optics 
lj. Processes of that type known so far can be sepa- 
rated into two groups. The first group comprises differ- 
ent interference effects, such as two-beam interference 
of photons in a beam-splitter In some schemes, the 
use of ancillary photons and conditional detection en- 
ables the realization of effective nonlinear interaction by 
means of linear optics Q . Besides their fundamental im- 
portance, such processes have tremendous potential for a 
quantum computing. In particular, such transformation 
of photon statistics as nonlinear sign shift 0,0] and arbi- 
trary photon number state filtering 4] have been realized. 
The manifestation of phase-dependent photon statistics 
in the mixed field of a narrow band two-photon source 
and coherent field has been reported in Ref. Q . Photon 
antibunching as well as photon bunching have been ob- 
served. One can expect that this scheme should be able 
to produce effectively sub-Poissonian photon statistics. 

The second group is formed by the processes in systems 
with nonlinearity and gain, where the transformation oc- 
curs due to the strong light-matter coupling pj. A ma- 
jor development is that of single-photon sources |(J. The 
strong light-matter coupling enables different methods of 
generation of the sub-Poissonian photon statistics, such 
as micromaser Q , single-atom resonance fluorescence , 
single-trapped-atom laser , and cavity-QED laser [!j . 
In Ref. , a method has been proposed of the transfor- 
mation of a vacuum state of cavity mode into arbitrary 
quantum state of light by means of succession of strong 
atom-field interaction processes inside the cavity. 

In the present paper, we present an opportunity of a 
new type: transformation of the photon statistics in a lin- 
ear homogeneous medium those refractive index depends 
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on the photon number distribution. As a physical realiza- 
tion of a medium with the photon statistics-dependent 
refractive index, we propose a collection of noninteracting 
quantum dots (QDs) embedded in a host semiconductor. 
This composite is called an excitonic composite |lllll2| . 

The optical and electronic properties of QDs are cur- 
rently in vogue, due to promise in semiconductor de- 
vice physics |l3j and, in particular, for the quantum 
information storage and processing 0, 0|. The ap- 
plication of QD s as potential quantum-light emitters 

EHEEGllii 

is now intensively discussed. In par- 
ticular, Ref. |19| demonstrates a quantum light source of 
multicolor photons with tunable statistics. The coupling 
of two QDs has been demonstrated as an effective means 
of emitted photons antibunching 20] . 

The peculiar property of a QD exposed to electro- 
magnetic field is the pronounced role of the dipole- 
dipole electron-hole interaction. Phenomenologically, 
this interaction can be introduced trough the local fields 
pH I22I I23L |24| . Obviously, in quantum electrodynam- 
ics the dipole-dipole interaction is due to the exchange 
between electrons and holes by virtual vacuum photons 
[2l| . Refs. 0, predict a fine structure of the absorp- 
tion (emission) line in a QD interacting with quantum 
light. Instead of a single line at the exciton transition 
frequency loq, a doublet appears with one component 
blue (red) shifted by a value Aui. The fine structure 
has no analog in classical electrodynamics. The value 
of the shift depends only on the QD shape, while the 
intensities of components are completely determined by 
the light statistics. In the limiting cases of classical light 
and single Fock states, the doublet reduces to a singlet 
shifted in the former case and unshifted in the latter one. 
In heterogeneous mediums comprising regular or irregu- 
lar arrays of QDs embedded in a semiconducting host, 
this mechanism is responsible for the photon statistics 
dispersion mentioned earlier. 

The electromagnetic response properties of a compos- 
ite medium comprising electrically small (of dimension 
much smaller than the wavelength) inclusions exposed 
to classical light can be modeled within the effective- 
medium approach modified 0, ^3 to include specific 
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properties of excitons as resonant states with discrete 
energy spectrum. To solve the problem of the quantum 
light interaction with a composite medium, an adequate 
quantization technique must be developed. Convention- 
ally, particular models of the classical light interaction 
with homogeneous media are exploited for this purpose. 
Nonhomogeneites are implemented into models by means 
of homogenization procedure on the classical stage. How- 
ever, the subsequent replacement of classical fields by cor- 
responding field operators often leads to a fundamental 
problem, which is the lack of correct commutation rela- 
tions for the field operators. To overcome the problem, 
some modifications of the quantization scheme have been 
proposed, such as, for instance, the noise current concept 
[2fi| . Nevertheless, in all cases to our knowledge, the con- 
stitutive parameters in classical and quantum optics are 
assumed to be identical. 

Here we propose an alternative approach: first we con- 
sider the interaction of quantum light with an isolated 
scatterer (QD) and then we develop a corresponding ho- 
mogenization tec hniq ue. The approach exploits the well- 
known technique [271 128| linking the macroscopic refrac- 
tive index to the forward scattering amplitude /(0) by 
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2 /(0) , 



(1) 



where k — lu/c is the vacuum wavenumber, c is vacuum 
speed of light, and p is the density of scatterers in the 
medium. Within the approach, the quantum field opera- 
tors satisfy the correct commutation relations automat- 
ically. Note that the approach is applicable not only to 
photons but to quantum fields of other origination (atoms 
and atomic nuclei, neutrons, etc.) and is in the basis of 
nuclear optics of polarized media [2j| . 

The paper is arranged as follows. In Sec. II we present 
a model of the single QD-quantum light interaction and 
evaluate the forward scattering amplitude operator. The 
fundamentals of quantum optics of excitonic composites 
— constitutive relations for field operators, planewave so- 
lution of Maxwell equations, input -output relations for 
a planar layer — are formulated in Sec. III. Sec. IV is 
devoted to the photon statistics dispersion. A general 
formulation of the transformation law for the light den- 
sity operator and the second-order coherence correlation 
function are set up and some particular examples of dif- 
ferent initial quantum states are considered. A summary 
of the work and outlook are given in Sec. V. The Ap- 
pendix provides some details of derivations presented in 
Sec. II. 



II. INTERACTION OF QUANTUM LIGHT 
WITH A SINGLE QD 

A. General relations 

We start with the evaluation of the forward scattering 
amplitude /(0) of quantum light interacting with a single 



QD. As has been pointed out by Cho 30], in the formu- 
lation of Hamiltonian of the interacting system " electro- 
magnetic field + condensed matter", the decomposition 
of the system into matter and field is ambiguous. In the 
approach called scheme [A] [3(J , only the transverse part 
of the total field — radiation field — is attributed to 
electromagnetic field. Charge carriers are assumed to be 
interacting via longitudinal near field which thereby is 
attributed to the continuum. In such a scheme, the free 
electromagnetic field is described in terms of the photon 
creation/annihilation operators fit, fi while the QD con- 
tinuum is described by the quantum field operators ^(r) 
and \&t(r). These quantum fields are fermionic [3(| and, 
consequently, satisfy the anticommutation relations: 



[*t( r ),$( r ')] + = 8(v-r>), 
[$(r),$(r')]+ - [§t(r),§t( r ')] + = . 



(2) 



Accordingly [3(|, the Hamiltonian of a QD exposed to 
quantum electromagnetic field can be written as 



H = "Ho + H ph + H-i + Ti-c ■ 



(3) 



Here Ho is the Hamiltonian of the charge carriers free mo- 
tion. The term " free" implies the carriers to be spatially 
confined in the QD but free of the influence of electromag- 
netic field and interparticle interaction. The Hamiltonian 



H 



hw fi^fi^ duj 



(4) 



describes photons in free space in the single-mode 
approximation The index ui ascribes the cre- 

ation/annihilation operators to photon of a frequency to. 
The creation/annihilation operators satisfy the commu- 
tation relation [S^, fi^,] = 6(oj—lu'). The light -matter in- 
teraction Hamiltonian 7i l in the length gauge Q] is given 
by 



Hj = -- I (V£o+£ a V)d 3 r, (5) 



where 



V(r) = e$ t (r)r$(r) 



(6) 



is the time-domain polarization operator, e is the elec- 
tron charge, V is the QD volume and £q is the free elec- 
tric field operator. Besides the Heisenberg representa- 
tion of operators £o(r) and "P(r), further we shell use 
their Fourier transforms with respect to time. In partic- 
ular, the electric field operator can be represented by the 
superposition £q(t) — £^ + \r) + £ ~\r), where super- 
scripts (±) mark contributions from positive and negative 
frequencies, respectively. Thus, the time-domain opera- 
tor £ q (r) is associated with the frequency-domain field 
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by the positive semi-axis Fourier transform: 

/>oo 

^ +) (r)= / Eo(r,w)dw. (7) 

J o 

The operators £q ^ (r) and £^ (r) are hermitian conju- 
gates. In the single-mode approximation 

Ihk 

E (r, w) = iy — e y a w e lkz , (8) 

where ej, is the unit polarization vector of free photons 
and A is the normalization cross-section in the xy plane. 
The Fourier transforms of the polarizability 'P(r) and the 
electric displacement f(r) are introduced in the same 
manner. 

The Hamiltonian 7i c in Q describes the interaction 
of particles in the QD. The retardation in that inter- 
action can be ignored in view of the electrical small- 
ness of the QD. Then, in the Hartree-Fock-Bogoliubov 
approximation, the Hamiltonian TL C is represented by 
H, c = H. C1 + 7i C2 , where 

H cl =e*J J $t (r) $ (r)( $t (r0 $ (r0> |^!£:, (9) 

v v 

K C2 =-z 2 JJ *V) *00 (*t(r)$(r')>^^ . (10) 
v v 

The term 7i ci accounts for the direct interaction of par- 
ticles while the term Ti. C2 defines the exchange interac- 



tion. Here and hereafter the angular brackets (...) stand 
for the operator's mean value. Note that terms @ and 
ljl comprise both static Coulomb interactions and local 
fields (dynamic Coulomb forces). As the next step, we 
reformulate the Hamiltonian Q in order to isolate the 
contribution of the local field. 

The quantum field operators \I/(r) and <3>^(r) are re- 
lated to real particles in QDs — atoms in crystalline lat- 
tice. The introduction of quasi-particles supposes aver- 
aging over the unit cell and transition to envelopes. In 
accordance with the envelope technique [3l| , the field op- 
erators can be represented by 

$(r) = y/v/N Y, in{R q )w ng (v) , (11) 

n,q 

where N is a total number of atoms in the QD, w nq (r) is 
the Wannier function in the q-th point of the QD crys- 
talline lattice, and R 9 is the radius-vector of that point, 
the index n stands for the collection of quantum numbers 
characterizing electronic states in the QD, and £ n (R g ) is 
the quantum field envelope satisfying fermionic anticom- 
mutative relations IjAip . 

Next we substitute Ijlip into and (|10p and carry 
out the averaging over the unit cell (see Appendix for 
details). Transition to the continuous field of envelopes is 
realized by means of the replacement (|A1 1|> : as a result, 
the Hamiltonian Ti c is expressed by the superposition 
H c = H oc + AH, where 



H 



nc 



e 'E // k>K>) (L(r')L(r')) -&(r)£n(r') (&(r)£„(r')) 



d 3 r d 3 i 



r - r' 



(12) 



vv 



AH = o E // /wG(r-r> w [|t(r)Mr)<^^ (13) 



m,m y y 



Here G(r) is the electrostatic Green's tensor of free space 
and li nn ' is the transition dipole moment (see Eq. (|A9p 
and following comment). The term 112[) describes the 
influence of static Coulomb forces while the Hamiltonian 
A7Y (|13p accounts for the local fields. 



B. Model Hamiltonian 

For further analytical progress, we use a set of simpli- 
fications. First, we restrict the consideration to two level 
model of exciton in QD, with \g) and |e) as its ground 
and excited states, and neglect intraband transitions, i.e., 



we let ^inn = 0. Exciton is assumed to be strongly con- 
fined in the QD; under that assumption the wavefunc- 
tions of states \g) and |e) are the same ;.3!j- Denoting 
them as £(r) we present the quantum field envelopes by 
£g,c = a g ,c£( r ) and e = a| e £*(r), where a| e and a gje 
are the operators of the electron creation and annihilation 
in the ground and excited states, respectively. Then we 
omit the term 7i c which accounts for the Coulomb in- 
teractions and is small in the strong confinement regime 
|3l|. Lastly, in order to eliminate fast oscillations, we 
make use of the rotating-wave approximation ,1] . As a 
result, we come to the light-matter Hamiltonian © with 
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H c = AH and 



(14) 



H Y = H J '(g u a + a u + 9^a_alj du , (15) 



AH = hAuj{a_{a + ) + , 



(16) 



where loq — (e e — e s )/h is the resonant frequency of the 
transition in QD and e giC are the energy eigenvalues in 
the ground and excited states, respectively, a z ± are the 
Pauli pseudospin operators. The quantity 



9u 



-i(n-e y ] 



hcA 



(17) 



is the coupling factor for photons and carriers in QD; 
(j, = (/x e g ' e y) e y = (A*g e ' e y) e y The depolarization shift 



47T 

Alu = —(fi-Nfj.) 



(18) 



characterizes the role of local fields inside QD. For a QD 
exposed to a classical electromagnetic field, the depolar- 
ization shift Auj has been introduced in Refs. [32. l33|. 
The depolarization tensor is given by 



V 



- v =--^ j j ie(r)ne(r')rv r ®v r 

V V 



1 



r - r' 



7T H'rdV, 
(19) 

where V r <g> V r is the operator dyadic acting on variables 
r. Physical meaning of the local-field correction becomes 
apparent if we rewrite the total interaction Hamiltonian 
Hi + Hq = Hi + AH in the form as follows: 



Hj + AH = — - j {Ve L + £ L V) d*r , (20) 



where the local field operator £l is related to the free 
electric field operator by 



H.c. 



(21) 



This relation demonstrates that the local field manifests 
itself as the QD depolarization, i.e., the screening of ex- 
ternal field by charges induced on the QD surface. The 
depolarization field defined by the second term in the 
right side of Eq. (|21|l is a classical field and is expressed 
by a C number. The result follows from that the depo- 
larization field in the quasi-static approximation is lon- 
gitudinal and thus is not quantized. That is, this field is 
not expressed in terms of photonic operators and, con- 
sequently, does not experience quantum fluctuations. A 
similar situation has been considered in Ref. |84j | under 
the analysis of Lorenz-Lorentz correction for quantum 
light in optically dense media. 



Any QD is essentially a multilevel system. However, 
the joint contribution of all transitions lying far away 
from a given resonance can be approximated by a non- 
resonant relative permittivity Sh- The host semiconduc- 
tor's relative permittivity is also assumed to be equal 
to eh- For analytical tractability, let £h be frequency- 
independent and real-valued. This allows us to put 
£/j = 1 without loss of generality. The substitutions 
c — > cj \fehi I 1 ~~ > AV 1 y/£h m the final expressions will 
restore the case ^ 1. 



C. Equations of motion 

In the previous sections, we have formulated the light- 
matter interaction Hamiltonian © and defined all its 
components in terms of the operators a u> aj, , &± and a z 
by Eqs. I|14 |l — Qlfl jl . In the Heisenberg representation, 
the equations of motion for operators a w , a _ and <x z are 
given by 



d 



iuJo&_+a z / g^a^ dcu + iAiv(j z (&_) , (22) 



(23) 



d . 

— a w = -lua^ - ig^cr^ , 
at 

oo 

d f 

—a z = 2 g u (ala_ - <r + a w )dw 



+2fiAw(d-_((T + ) + &,(&_)) , 



(24) 

Let at the initial time t = — oo the QD be in the ground 
state \g) ((a z (—oo)) = —1) and be exposed to an ar- 
bitrary quantum state of light. The interaction is as- 
sumed to be weak, i.e., only first-order terms in g u are 
of importance. From Eq. 122|) follows that a± and 
(<7 ± ) ~ 0{gu), which means that da z /dt ~ 0(|<7w| 2 ) 
in 124fl . Consequently, the weak interaction model al- 
lows us to put a z (t) — a z (~oo) in the further anal- 
ysis. Since the final expressions will obligatorily in- 
volve averaging over the ground state, i.e., the quantity 
(a z (t)) w ((T z (— oo)) = —1, it is convenient to make the 
replacement a z — > — 1 straightaway in Eq. Q22[l. As a 
result we come to 

oo 

d f 

— (&_) = -i(u + Au))(a_) - i j g u (a u ) dui . (25) 

o 

The integration of this equation with respect to time 
yields 

oo t 

(a_) =-i Jg u J (a u (t')) e -^°+A")(*-* V du . (26) 

-oo 

The substitution of lt2*B|) into 1(2*2^) gives us 

oo 

d f - 

— &_ = -uv &_ - i I g^Quj did , (27) 
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where 



e u (t) = + / (a w (t')) e -«(^+A W )(i-t') dt ' . (28) 



By integration of 127(1 we obtain the explicit relation 

oo t 

o-_ = -t ^ /" 6 w (t') e-^-t'Ut' duj (29) 

-oo 

between the operators a u and <r_, which defines the in- 
teraction of the QD-exciton with quantum light in the 
weak-field limit. This equation serves as a basis for the 
further analysis because it allows us to develop the cor- 
relation between the electric field and the polarization 
operators in QD. Going to the frequency domain and uti- 
lizing the relation P(r) = |£(r)| 2 /x*<r_ + H.c. and Eqs. 
© , (|2*5|) and it2*9"|) we come to the frequency-domain po- 
larizability operator as 



P(r) 



\»\ 2 /hV 

UJ — UJq + iO 



E (r) 



Auj (E (r)) 
u> — ujo — Aw + iO 



(30) 

with Eo(r) given by Eq. JHJ). 

Let us restrict further consideration to monochromatic 
fields. This allows us to omit from here onwards the index 
u> in operators a w and aj, and omit explicit mention of the 
frequency dependence in complex-amplitude operators. 



D. Forward scattering amplitude 

Let the field Eo(r) be scattered by an isolated QD. 
Obviously, the total field E(r) outside the QD is a su- 
perposition of the incident and the scattered fields. In 
terms of the polarization induced in the QD, this field is 
expressed by the relation 

E(r) = E (r) + (VV • +fc 2 ) j P(r') y^d 3 r' , (31) 

v 

which follows from the Maxwell equations. In turn, the 
induced polarization is related to a weak incident field 
by the linear relation (|30|l . Substitution of (|30|l into 1(31(1 
allows to couple the incident and the total field, that is, 
to define the scattering operator |35T ]. 

The analogous relation with operators replaced by cor- 
responding C-numbers holds for classical fields. For elec- 
trically small scatterers, the coupling between the po- 
larization and the electric field is local. This allows the 
introduction of the forward scattering amplitude by 



P(r o ) = p/(0)E(r ), 



(32) 



where ro is the scatterer radius-vector. In the far zone, 
the scattered field is a diverging spherical wave; thus, 



E(r) = E (r) 



V 



(VV 



_ fc2) e^M /(0)Eo(ro) _ (33) 
l r _ r o| 



Then, the forward scattering amplitude / (0) is found to 
be proportional to the QD polarizability [21| and is given 
by 



/(0) = - 



i 



hV u> — lu — Aw + iO 



(34) 



Going to quantum light we cannot restrict ourselves to 
the formal replacement of electric field and polarization 
in ((32(1 and 1)33(1 by corresponding operators, keeping for- 
ward scattering amplitude in the form of (13411 . As follows 
from Eq. (|30() . because of the local field impact, the po- 
larization operator is represented by a superposition of 
the electric field operator and its mean value. As a re- 
sult, Eq. (|30|l can not be reduced to Eq. (|32|) written for 
operators. 

The necessary generalization of (|32(l exploits represen- 
tation of the forward amplitude as an operator in the 
space of quantum states of light: 



P( r o) - [^°) £ °( r °) + Eo(r )/(0) 



(35) 



As one can see, the product of the operators /(0) and Eq 
enters this equation in symmetrized form, analogously 
to the interaction Hamiltonian (JSJ. Using ©, l|3U|l and 
(|35|l . we derive the equation 

/(0)fi + a/(0) = 2/ 2 (0)a + 2(o)[/i(0) - / 2 (0)] , (36) 

for the forward scattering amplitude operator /(0), where 



/i, 2 (0) = - 



k 2 \»\ : 



1 



HV to — u>\ 2 + iO 



(37) 



are the partial forward scattering amplitudes with uj\ — 
luq + Aoj and <x>2 = ^o- To solve ()36|l with respect to /(0), 
we introduce the operator 



\2m+ l)(2m| 
' V2m + 1 



(38) 



where |n) stands for the n-th Fock state. The operator 
satisfies the identity 



£a + a£ = 1 



(39) 



We snow seek solution of (|36|) in the form of /(0) = 
X£ + Y. Substituting this representation into Eq. 1(36(1 
and equating the terms of the same order in a, we find 
coefficients X and Y and, thus, write the forward scat- 
tering amplitude operator explicitly: 



/(0) = 2<a)[/ 1 (0)-/ 2 (0)]| + / 2 (0). 



(40) 



The total electric field operator is expressed by 1(31(1 with 
the polarization operator defined by ((35(1 . 



6 



III. QUANTUM OPTICS OF EXCITONIC 
COMPOSITE 



A. Constitutive relations for field operators 

Consider a heterogeneous medium comprising a cubic 
lattice of identical spherical QDs embedded in a semicon- 
ducting host. The QD radius and distance between QDs 
are assumed to be much less than the wavelength in the 
host medium. The optical properties of such a composite 
medium are isotropic and can be described by an effective 
index of refraction. Let the optical density of the com- 
posite be small enough to neglect the electromagnetic 
coupling of QDs, i.e., to neglect the Mossotti-Clausius 
correction. Since the forward scattering amplitude turns 
out to be an operator, the optical response of the compos- 
ite is described, instead of (JTJ , by the index of refraction 



« 2 = 1 + ^/(0) 



(41) 



which is an operator in the space of quantum states of 
light. Obviously, the same is relevant to the relative 
permittivity e = n 2 . Utilizing the well-known relation 
D = E + 47rP and l|35|) , we come to the operator consti- 
tutive relation 



D = -(fE + Ee), 

with the effective operator permittivity 
£ = 2{a)in\ + (1 - 2(a)£)n 



2 

2 ' 



(42) 



(43) 



where n\ 2 = 1 + 47r j o/ 1 , 2 (0)/fc 2 . The constitutive re- 
lations (|42|l and (|43|) satisfy the principle of correspon- 
dence: in the classical limit D — > (D) we come to 
(D) = n.f(E), i.e., the correct constitutive relation for 
an excitonic composite exposed to classical light 0, . 
Correspondingly, in that case the forward scattering am- 
plitude is a C-number and is expressed by (j34(l . Even 
so, the local fields may result in physically observable 
effects, such as polarization-dependent splitting of the 
gain line [32l l33T | in asymmetric QDs. The splitting was 
observed experimentally in [3fl l37| and was referred to 
as a manifestation of the exchange interaction. Let us 
emphasize once more time that here we deal with the 
same effect presented in phenomenological and micro- 
scopic languages. 

Because of l(39f) , the operators e (|43|) and d do not com- 
mutate and, consequently, the symmetrized representa- 
tion of the constitutive relation (|42ll is of fundamental 
importance. Note that the operator e is non-Hermitian 
even in the absence of dissipation. Its physical meaning 
we discuss later on. Here we only stress that the non- 
hcrmitian nature does not result in any physical contra- 
diction since the operator e does not represent any ob- 
servable quantity. In contrast, the electric displacement 
operator T>, related to its Fourier amplitude D by the the 



procedure analogous to Eqs. J7J, is Hermitian as it must 
be for observable quantities. The hermicity is due to the 
symmetrized notation on the right-side part of Eq. H42fl . 



B. Generalized plane waves 

Consider the propagation in the excitonic composite 
of free (transverse) electromagnetic waves. We assume 
n 2 2 > 0; thus we exclude from the analysis narrow for- 
bidden zones between longitudinal and transverse reso- 
nances in the vicinity of frequencies Wi,2- The propaga- 
tion is described by the ID frequency-domain operator 
Maxwell equations 

ikH = e z x— , — _( e E + Ee)=e, x — . 44 
oz 2 oz 

This system can easily be transformed to a set of wave 
equations for the operator E — (E) and the mean value 
(E): 



d 2 



+ k 2 n 2 )(V) 



(^2+^i)(E-<E))=0. 



(45) 



(46) 



The elementary solution of these equations leads to the 
expansions 



/n 2 



(a -(a)) 



(47) 



H 



a ikni z 



ni(a) + 



Akri2 z 



n 2 (a 



(48) 

for the mode traveling in the direction z > 0. Eqs. I(47|) 
and H48JI display the local field effect in the excitonic com- 
posite: decomposition of a single mode into coherent and 
incoherent parts with different refractive indexes. Thus, 
as quantum light propagates through the excitonic com- 
posite, the light statistics is altered and spatial modula- 
tion of the amplitude with the period 



Ln = 



2tt 



ni - n 2 \k 



(49) 



occurs. The period is strongly dependent on frequency 
due to the frequency dependence of partial refractive in- 
dexes n\.2 determined by Eq. 1)37(1. In turn, the alter- 
ation leads to transformation of the light quantum coher- 
ence of the second and higher degrees. To some extent, 
the situation is similar to that in nuclear optics |29j | : the 
existence of two spin-dependent refractive indexes of the 
particle's wavefunction causes the spin rotation as the 
particle travels through a polarized medium. 

Physically, the result obtained can be understood as 
follows. There are two mechanisms of quantum light scat- 
tering by a QD. The first one - quasi-classical - is due 
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to the induction in the QD of the observable polariza- 
tion and, therefore, it persists in the classical limit. The 
corresponding frequency turns out to be shifted by the 
value Aw. This mechanism changes the QD state due 
to the local-field induced electron-hole correlations and 
thus provides inelastic channel of the scattering. The 
first term in the right side of l|47|l conforms to this chan- 
nel. The aforementioned non-hermicity of the operator 
e even in the absence of dissipation is due to inelastic 
scattering. The second mechanism manifests itself in the 
scattering of states of light having zero observable field 
(e.g., a single Fock state). Such states do not induce 
observable polarization and, consequently, any frequency 
shift. This mechanism has no classical analog and pro- 
vides elastic channel of the scattering. The presence of 
the inelastic scattering channel signifies the change of the 
scattered light's quantum state, i.e., the photon statistics 
dispersion. Thus, the photon statistics dispersion and the 
resulting phenomenon of the electromagnetic field spatial 
modulation are the direct manifestations of the fine struc- 
ture of excitonic line, inherent in an isolated QD |2lj . 

Going to time-domain representation of operators 147|) 
and igSJl, £(r, t) and ?i(r, t), by the procedure defined by 
Eq. (0, we can show that the condition 



J (£(r, t) x ?i(r, t))e z dx dy — const 



(50) 



holds true for an arbitrary quantum state. Physically, 
this condition reflects the momentum conservation for 
fields determined by operators (|47|l and l|48|l : a mean 
value of the energy flux through an arbitrary closed sur- 
face is equal to zero for an arbitrary quantum state of 
light. That is, dissipation is absent in the medium in 
spite of the non-hermicity of the relative permittivity op- 
erator e. The conservation low l|50|l testifies to physcal 
correctness of electromagnetic waves defined by the op- 
erators and (|4"g|) . 



C. Input— output relations 

Consider a planar layer of the QD-based composite ma- 
terial of thickness d located at \z\ < d/2 exposed to a 
quantum field. The problem to be solved is to ascertain 
relations between operators of incoming, reflected and 
transmitted waves. For that, we make use of Eqs. lB7|) 
and H48f) and utilize the layer's reflection R and trans- 
mission T amplitudes (see, e.g., Eqs. (6.130)-(6.132) in 
Ref. |26]). The reflection and transmission amplitudes 
satisfy the property |i?| 2 + |T| 2 = 1 which is a reflection 
of the scattering matrix's unitarity. Obviously, the total 
field outside the layer can be represented by 26] 



e = e 1 : 



E° 



E 11 



E° 



(51) 



where E'_£ are the operators of incoming waves propa- 
gating along the positive (index "+") and negative (in- 
dex "— ") directions of the axis z. The superscript "out" 



marks outgoing waves. By analogy with (JSJ) the incoming 
waves are given by 



E" 




A::tkz 



(52) 



where creation/annigilation operators a±/a± satisfy or- 
dinary bosonic commutation relations. 

In view of the commutativity of operators On the 
right sides of Eqs. iHTj) and the layer's transmis- 

sion/reflection balance can be written separately for each 
of waves exp(ikni, 2 z): 



Ti(a±) + i?i(a=, 



(53) 



-out 

a ± 



(O = T 2 (a± - (o±» + R 2 (a T - (a T )) . (54) 



Here T\ 2 and are the transmission and reflection 
coefficients of layers of the thickness d and refractive in- 
dexes n\ and n 2 , respectively. Using 15311 154|) . we can 
represent outgoing fields as 



e^i 



-Likz 



[T 2 d± + R 2 a T 



+ (T 2 - Ti)(a±) + {Ri - R 2 )(d T )} . (55) 

It can easily be shown that, in view of |i?i,2| 2 + IT1.2I 2 = 
1, the outgoing electromagnetic field satisfies the correct 
commutation relations for field operators and the energy- 
momentum balance condition for incoming and outgoing 
waves for light in arbitrary quantum state. 

Equations (|55|l show that the quantum light transmis- 
sion through a layer of excitonic composite leads to the al- 
teration of light coherence. For example, using the obvi- 
ous identity ((a±-(a±)) t (a±-(a±))) = (a±a±)-|(a±)| 2 , 
we derive 



(|E- 



out|2\ 



I To 



ITJ 2 



|T 2 | 2 )|(a + )| 2 . (56) 



The photocurrent measured by a detector located in the 
region z > d/2 is proportional to (|E° ut | 2 ). Thus, the 
photodetector signal is no longer completely determined 
by the mean number of photons (first term in (|56|l 'l but 
depends also on the light statistics (second term in J5BJ)). 
This effect can be exploited for experimental verification 
of theoretical predictions of the present article. The dis- 
tinction of the reflection/transmission properties of the 
excitonic composite layer for coherent and non-coherent 
components can be used for the design of new quantum - 
optical light transformers. 



IV. PHOTONIC-STATE DISPERSION 

The complete characterization of a quantum state of 
light is given by the density operator. In this section, we 
evaluate this operator in an excitonic composite medium. 
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Along with that, the issues of interest are the light 
characteristics measurable in typical photodetection pro- 
cesses. As an example of such characteristics, further 
we consider the normally-ordered normalized time-zero 
second-order correlation function g( 2 '(z) of the excitonic 
composite, which characterizes the quantum-mechanical 
second-order degree of coherence at a given cross-section 



A. Density operator 



Let us rewrite Eq. (|4T|) as 



where a(z) = a + k[z) and 



k(z) 



n 2 i(n 1 -n 2 )kz 



(57) 



(58) 



The quantity k(z) is a periodic function of the spatial 
coordinate z with oscillation period L 149(1. Operators 
a(z) and a*(z) satisfy the ordinary bosonic commutation 
relations. That is, they can be treated as annihilation 
and creation operators of the quantums of some bosonic 
field. Utilizing then operator identities (3.47) and (3.48) 
from 261, we come to 



a = D\aD K , 



where 



D K = exp(Ka^ — n*a) 



(59) 
(60) 



is the displacement operator with k as the displacement 
factor. Hereafter we omit the ^-dependence. According 
to (j59*|) , as light goes through the excitonic composite the 
spatial evolution of the light density operator is governed 
by the unitary transform 



DlgD K 



(61) 



where g is the initial operator of the light density induced 
by a source at the input. 

The distinguishing property of the transform (|61|l is 
that it actually is not linear: the displacement factor k 
depends on the state of initial quantum light through the 
relation 



n=0 



+ l,n 



(62) 



Physically, this means that the field acting on the 
QD-exciton is self-consistent in the Hartree-Fock- 
Bogoliubov sense (see Eqs. (fl^|) and As a result of 

self-consistency, nonlinearity of the quantum-mechanical 
motion of carriers in QD arises; see the second term on 
the right side of Eq. 



The operator transformation law (|61|l can be rewritten 
in terms of familiar matrix elements. As the first step we 
represent g' by the expansion 



?' = 2J QmnD- K \m){n\DL K . 



(63) 



In accordance with Ref . |26j , 



D-Jm) = exp 



Y,S mp (n)\p), (64) 



where 



*S'mp( K ) — 



{ - K) P-m l^_ L (v-™ )(| K |2 )j m<p ^ 



(65) 

and Lm\x) are the generalized Laguerre polynomials. 
Then, using H63fl - I|65[l we can proceed from i|61|) to the 
law of transformation of matrix elements: 



CXpf - \ y ' grnnSmp{ K )S nq (n) , 



(66) 



Equations (|61|l and (|66(l determine the transformation 
law for quantum light with arbitrary statistics. Consider 
several particular cases as examples. For a single Fock 
state \n), (a) =0 and k = 0. Taking Dq = 1, into ac- 
count we conclude that the excitonic composite does not 
transform single Fock states. Analogous situation oc- 
curs for thermal light where g mn ~ S mn and, accordingly 
(|62|l . (a) — 0. The spatial evolution of a coherent state 
1/3) is characterized by the transform g' = D\\f3) ((3\D K = 
D- k D-p\0)(0\d' < _ (} D^ k . Taking into account the well- 
known identity |2(j D- K D-p = D-p_ K exp[Tm(K/3*)], 
we come to 



\k + 0)(k+0\ 
— 3 e*'™ 1 ~™ 2 )' cz 



— 3 e i ( ni ~~ ri2 ) kz 



(67) 



This result implies that initial coherent state is trans- 
formed into another coherent state with modified coher- 
ent amplitude. The coherent amplitude of transformed 
state experiences spatial beating as light passes the com- 
posite (factor exp[i(ni — n2)kz] in l|67l) ). The mean pho- 
ton number in the composite y/ ri2/ni|/3| 2 holds invari- 
able and significantly differs from that the source pro- 
duces in vacuum, \3\ 2 . The alteration is determined by 
the ratio ni jn\ which varies in a broad range depending 
on the photon energy. There are different energy bands 
in the excitonic composite, where the classical source ef- 
ficiency either rises (ni > n-y) or falls (n-i < ni), in com- 
parison with the source efficiency in vacuum. 

As one more example, consider the transformation of 
the Fock qubit. Let |*) = /? |0) + 0i\l) , where the co- 
efficients /3q,i are coupled by the ordinary normalization 
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| A) | 2 + |/3i| 2 = 1 • According to © and the trans- 
formation is given by 

q> = \V>){y'\ = D- K \*){*\Dl K , (68) 

with (a) = /3q/3i m lHH f° r the displacement fac- 
tor k. Obviously, in the coherent-state basis \^>') = 
(PoD- K + PiD-Ka^lO) . Taking then into account the 
relation [tf,D- K ] = -dD- K /da = -k*D_ k (see [H Eq. 
(C.17)]), we obtain |*') = (fig + x k* + ^)\k) '. Then, 
proceeding to the Fock-state representation we come to 
the expression 

I 12 00 n 

W) = exp(-^-)£ A_ + filK* + \n)- (69) 

Equations (|68|) and l|69|) define the spatial evolution of the 
Fock qubit in excitonic composite. Let us stress that l|6^|) 
includes a full set of Fock states. An essential physical 



conclusion can be deduced from that: new states absent 
in light in the beginning arise as light travels through the 
excitonic composite. 



B. Second— order coherence 

To illustrate the effect of the photon statistics disper- 
sion, we calculate the second-order correlation function 
g( 2 \z), which is represented in terms of the field opera- 
tors by 

g^ 2 \z) = (E^E^EE) / (E^E) 2 . 

Let us consider light prepared in an arbitrary quantum 
state |$). Then, after some algebraic manipulations, we 
obtain the second-order correlation function for the field 
operator ijiTjl 



(2) J^n W{z)A n _ 1 + 2k(z)V7L4» + y/ n(n + l)A,+i| 2 



E„^K| 2 +(^-l)|(a)| 2 



1 2 



(70) 



where A n — (n\^). Without taking local 
fields into account, the correlation function 
g^(z) coincides with the correlation function 

$ ] = £„^ + i)K+il7(£„™KI 2 ) 2 of the 

incident field Eo. The distinguishing property of 
g( 2 '(z) is its frequency dependence which is due to the 
pronounced frequency dispersion m t 2 — ^1.2 (w) of the 
excitonic composite. Also, the second-order correlation 
function experiences the spatial modulation with period 
Lo. 

The most intriguing peculiarity in the photon statis- 
tics appears in the vicinity of the resonance uj = W2- For 
the light states with a large enough (a), the inequality 
n 2|(a)| 2 1 holds true because of the growth of n-i in 
the vicinity of the resonance. As follows from (|70l) . in 
that case the spatial modulation becomes a high-order 
infinitesimal and g^ 2 \z) — > 1 independently of the light 
state 1^); that is, the photon statistics tends to Poisso- 
nian (completely uncorrelated photons) as light propa- 
gates through the excitonic composite. This means that 

the excitonic composite reduces the photon correlation 

( 2) 

degree for photons with both positive (<?q > 1) and 
negative (g^ < 1) correlation. Indeed, g^ 2 '(z) < g ' in 

the former case, whereas g^ 2 \z) > g^ in the last one. 
Another situation arises for a pure Fock state |n) with 
(a) = 0. In this case Eq. (JZD) gives g^ = l-l/n. That 
is, the correlation function l|7U|l coincides with the free- 



space correlation function. In general, negligible chang- 
ing of the correlation function takes place for all quasi- 
Fock states characterized by the condition n2\(a)\ 2 <C 1. 

A spatial modulation of the second-order correla- 
tion function of light is characteristic for the Hanbury- 
Brown-Twiss (HBT) interference Q of quantum light 
emitted by two sources. As different from the HBT 
effect, the spatial modulation g^ 2 '(z) H7U|) in the excitonic 
composite arises as a result of the interference of coher- 
ent and incoherent components of the field. The trans- 
formation of the statistics of quantum light interacting 
with a finite-thickness dielectric slab with a C— number 
relative permittivity was considered in Ref. 26]. The 
statistics transformation is provided by the light reflec- 
tion/refraction at boundaries, similar to that as a disper- 
sionless dielectric slab modifies the frequency spectrum 
of a signal due to the interference of waves reflected from 
boundaries. Another effect of such a type has been pre- 
dicted in Ref. 38j: appearance of spatial quantum cor- 
relation in light propagating through a multiple scatter- 
ing random medium. The quantum correlations depend 
on the state of quantum light and are due to the cross- 
correlation between different output modes. In our case, 
the physical mechanism of the transformation is com- 
pletely different: we predict the transformation effect to 
be inherent to a homogeneous medium. That allows us 
to identify the effect as a new type of optical dispersion. 

It should be noted that the notions " linear regime" and 
"weak coupling regime" as applied to excitonic compos- 
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ite mediums are not identical. In the presence of local 
field, the polarization 130|) is linear in the incident field 
(linear regime) but contains a quadratic term in the os- 
cillator strength, i.e., the term 0(|/x| 4 ). Nonlincarity of 
that type violates the weak coupling regime and leads 
to the transformation of the photon statistics similar to 
that takes place in the standard Jaynes-Gummings dy- 
namics 1]. In other words, in the presence of local field 
the light-QD interaction is characterized by two coupling 
parameters: the standard Rabi frequency and the new 
depolarization shift Alu. 



V. SUMMARY AND OUTLOOK 

In the present paper we have presented a theory of 
the quantum light interaction with excitonic composites 
whose relative permittivity is predicted to be an oper- 
ator in the space of quantum states of light. As a re- 
sult, the medium's refractive index turns out to be de- 
pendent on the photon number distribution providing, 
by analogy with the frequency dispersion and the spa- 
tial dispersion, the photon statistics dispersion of light. 
Self-organized lattices of ordered QD-molecules and ID- 
ordered (In,Ga)As QD-arrays [3!J,|4(| serve as examples 
of high-quality excitonic composites. 

Light transformation in the excitonic composite is car- 
ried out by means of the displacement operator — one 
of the canonical transforms of quantum optics [l| . How- 
ever, as different from ordinary cases, the displacement 
factor k in an excitonic composite depends not only on 
the photon energy and momentum (in mediums with spa- 
tial dispersion) but also on the photon statistics. In par- 
ticular, an excitonic composite provides different prop- 
agation regimes for coherent and non-coherent compo- 
nents of total field. That is, the frequency response char- 
acteristics of excitonic composite-based structures may 
differ in coherent and non-coherent light. For example, 
a distributed Bragg reflector with layered components 
made of excitonic composite materials would manifest 
shifted reflection bands for coherent and non-coherent 
light. In the same manner, an excitonic composite-based 
microcavity will possess, depending on the light statis- 
tics, two alternate sets of eigenfrequencies. An interfer- 
ometer constituted by two collinear plane mirrors of exci- 
tonic composite material can serve as the simplest model. 
The excitonic composite provide a new (not related to 
anisotropy) birefringence mechanism: coherent and non- 
coherent components of a plane wave obliquely incident 
on plane excitonic composite interface will propagate at 
different angles. Equations (|53l) and (|54|l shows that an 
excitonic composite layer can serve as a beam splitter 
with a modified law of transformation of the field oper- 
ators. Along with the second-order coherence transfor- 
mation, the optical effects mentioned create a potential- 
ity for experimental verification of the predicted photon- 
statistics dispersion. 

The theory elaborated in the present paper deals with 



the single-mode case that implies the single pair of 
creation/annihilation operators sufficient to characterize 
the fields Q47|l and (|48|l . In the multi-mode case each 
mode will be characterized by its own mode-number- 
dependent value (a) and, correspondingly, permittivity 
operator. As a result, then the photon statistics dis- 
persion will be accompanied by spatial dispersion. In 
classical crystal optics [41| , spatial dispersion necessarily 
entails optical anisotropy even in structurally isotropic 
mediums. The preferential direction is determined by the 
propagating wavevector and refractive indexes are dif- 
ferent for transverse and longitudinal field components. 
In an isotropic medium with the multi-mode photon- 
statistics dispersion we meet another situation: in spite of 
the spatial nonlocality the permittivity operator is scalar. 
The situation can be understood from the following. A 
preferential direction is absent in an excitonic compos- 
ite since it can be related only to a physically observable 
wave. In contrast, in the case considered the spatial dis- 
persion mechanism concerns with operators and disap- 
pears on proceeding to observable fields. 

The excitonic composite considered in this paper is 
an example that demonstrates the existence of mediums 
characterized by the operator constitutive relation i|42|) 
and, consequently, displaying the photon statistics dis- 
persion (|43|l . It would be of interest to search other artifi- 
cial (and, possibly, natural) materials with photon statis- 
tics dispersion and the statistics transformation law not 
necessarily defined by il4"5|) and respectively. Since 
the photonic state dispersion in excitonic composites is 
due to dipole-dipole interactions in QDs, one can ex- 
pect the dispersion of that type in mediums with strong 
dipole-dipole interaction, such as Bose-Einstein conden- 
sates |25| and ultracold atomic ensembles |42j| . 

The developed model of quantum-optical properties 
of excitonic composites leads to a set of functional equa- 
tions, in which both quantum fields to be found and pre- 
determined coefficients are operators in the same space 
of quantum states. The investigation of general prop- 
erties and methods of solution of such equations con- 
stitutes an intriguing mathematical problem. One can 
expect that the equations of that type will find applica- 
tion in studying of quantum fields of different physical 
origin. The matter of special interest would be the case 
when the equation solution and the coefficients are non- 
commutative — as it takes place in our paper. 
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APPENDIX A: ENVELOPE METHOD FOR THE 
LOCAL FIELD-INDUCED CONTRIBUTION TO 
THE HAMILTONIAN 

The application of the envelope representation (|11|) to 
the Hamiltonian He given by Eqs. 10 and (|10l) makes 
possible averaging of © and IjlUII over unit cell and pass- 
ing to continuous field of envelopes, i.e., allows to rewrite 
He in terms of 1|12H and l|13fl ■ In spite of the fact that the 
procedure is well known (see, e.g., |43|\ its realization in 
our case contains some peculiarities and thus requires a 



more detailed presentation, which is given for the com- 
ponent Hci in the following. The case of Hc2 can be 
considered by analogy. 

First, we stress that Eqs. (0 and (|1 l|l produce the 
commutation relations for envelopes: 

(Al) 

[£ l (R,),^(R 9 0]+ = [&(R 9 ),&(R 9 <)]+ = 0. 

Substitution of (fTTJl into leads to 



H C1 = ( ^ ) E E^( R P)^'( R P') (ii(^AK q '))K P n't m ' , (A2) 



N 



m,m p,p 
n,n' q,q 



where 



K vvm - N 2 I I < P ( r )"W ^W mq {v')w m , q , (r') 



nn'mm' 

v P v„ 



and V p is the volume of the p-th unit cell. 

For evaluation of integral (| A3|) we introduce vectors 



x = Hp — H q and y = r — r' and take into account that 
y <C x within the unit cell. This allows us to represent 
the denominator of (|A3(I by the expansion 



y v, 



In the standard envelope technique, the third term is sup- 
posed to be a small correction to the Coulomb interaction 
and conventionally neglected |43j |. But, in our paper we 
demonstrate that the correction can be of importance 
and may result in qualitatively new effects. 

Substitution of (|A3|) into i|A4|l leads us, in view of the 
orthonormality of the Wannier functions w np (r), to 

T^pp'qq' S nn iS mm iSpp>Sqqi N 

^nn'rnm' ~ |T> tTI ' TT^nn'pp' 1 



V X ®V X ' 



1 

Rp — R<j I 



G 



mm' qq' , 



where 



G„ n 'p P < = y rw; p {v)w nlp , (r) d 3 r . 



(A5) 



(A6) 



For further consideration it is convenient to turn from 
Wannier to Bloch functions by means of the procedure 



(r) = — L E exp(ikR p )u„ fc (r) . (A7) 



Using this representation we reduce Eq. (|A6|) to 

Gnn'pp' = —Tf E ex P[*( k/R p' ~ kRp)]/x„ n / fcfc - , (A8) 



where 



Mi 



'kk' 



eN 



<(r) d 3 



(A9) 



is the matrix element of the dipole transition in a 3D 
medium the same as the QD material. As a rule, the most 
important contribution into dipole moment is from direct 
transitions in the vicinity of the bottom of band |43j . 
For such transitions we can neglect the fc, k' dependences 
in the last equation and take fJ- nn 'kk' — l^nn' thereafter. 
Then, taking the identity J2k exp[ik(R p / — R p )] = NS PP >, 



(A10) 



we see that Eq. (|A6|) can easily be reduced to 



G 



nn'pp 



The concluding step is to turn to the continuous field of 
envelopes. For that we substitute l|A5|) and IjAlOfl into 

r and R„ — * r' 



Sll and perform the replacement R p < 
to get 



V v 
TV 2^ P 



d A r. 



(All) 



V 



G(r - r') . 



k.k' 
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